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2 ■ Abstract. The next generation mass probes will investigate DE nature by measuring non- 

■ linear power spectra at various z, and comparing them with high precision simulations. 

| Producing a complete set of them, taking into account baryon physics and for any DE state 

equation w(z), would really be numerically expensive. Regularities reducing such duty are 
essential. This paper presents further n-body tests of a relation we found, linking models with 
DE state parameter w(z) to const. -w models, and also tests the relation in hydro simulations. 

\ 1. Introduction 

' The background metrics of space-time can be written in the forms 

ds 2 = c 2 dt 2 _ a 2Q M 2 = fl2(r) (dr 2 _ 

defining ordinary or conformal time, t or r. In systems which abandoned cosmic expansion, the 
former expression yields the Minkowskian metrics: the no longer growing scale factor can be 
inglobated in spatial coordinates. In such systems t is the natural time coordinate; e.g., many 
cyclic motions become (quasi-)periodical in respect to t. On the contrary, when a system resent 
the overall expansion, the natural time coordinate is r. For instance, sonic waves in the pre- 
recombination plasma, are (quasi-)periodical in respect to r. This periodicity yields the setting 
of maxima and minima in CMB spectra, as well as the scales of BAO's in matter spectra. 

Let us consider then the comoving distance between today's observer and the LSB (last 
scattering band) , L = J rec di. It coincides with the conformal time T = f rec dr a photon takes to 
reach the observer from the LSB. Accordingly, if two models have the same comoving distance 
from the LSB, we expect them to share those features which are determined by the cosmological 
dynamics. As a matter of fact, it has been known since several years [T] that, if they share fi c j,, 
h and n s (density parameters, Hubble parameter in units of 100 (km/s)/Mpc, primeval scalar 
index), their linear spectra are quite close. In particular, models with the same conformal age, 
differing just because of DE state equation, exhibit almost coincident linear spectra. 



When non-linear effects are taken into account, the dynamics is no longer simply ruled by 
r. Non linear spectral deformations arise when matter halos form, abandoning the overall 
expansion. The point is whether they are so widely spread to yield a significantly discrepant 
spectral behavior at a given fc-scale. Greater k's correspond to scales entering a non-linear 
regime earlier; therefore, we expect greater discrepancies at greater fc's. 

A test of these expectations can only be performed in a numerical way, by running model 
simulations. Let us then outline that forthcoming tomographic lensing experiments (see, e.g., [2]) 
will be able to measure spectra with a precision 0(1 %) [3], so that this is the limit above which 
spectral discrepancies can no longer be disregarded. In turn, such experiments are expected to 
test spectra up to k ~ 10/iMpc _1 , corresponding to length scales 0(0.5 /t Mpc), well inside 
galaxy clusters. To test spectral regularities, then, hydro effects cannot be disregarded. 

A number of authors have considered these problems in some detail. The k value above 
which hydrodynamics causes spectral distorsions > 1 % has been found to be ~ 2-3 h Mpc -1 [3J. 
Models with an assigned polynomial DE state equation w(z) (models A) were then compared 
with models with constant w (models W, or auxiliary models), selected so to have the same 
L and the same a% (mean square density fluctuation at 8/i _1 Mpc), besides of Q& c and h [5]. 
Spectral discrepancies, within k ~ 2-3/iMpc -1 were found to keep within 1%. 

A critical feature of forthcoming mass probes is their capacity to explore the Universe at 
z > 0. Models with the same conformal age at z = 0, and different DE state equation, have 
different ages at z > and spectral discrepancies soon exceed 1 %. In principle, then, one should 
find an auxiliary model W(z) of A, at any z. But equal H at z ^ means different H at z = 0. 
Accordingly, simulations performed in boxes of identical size in units of h~ 1 Mpc, can no longer 
be compared (the suffix is used just here, to outline that h is taken at z = 0). If boxes of 
suitable different sizes are then used, sample variance is the dominant effect. 

This difficulty was overcome by abandoning the request that A and W(z) have the same 
H(z), and requiring them to share just the reduced density parameters w C) & = fl c ,bh 2 (weak 
requirement) [7J. The definition of p cr (critical density) however implies that, at any z, 

Q Cjb x H 2 = n Ctb x (8vrG/3)p cr , 

so that the uj c ^ tx £l c ^H 2 scales as Q c ,bPcr = Pc,b oc a -3 , indipendently of DE nature. Models 
with equal oj c ^ at z = will then share them at any z. The W(z) model expected to approach 
A at z, then depends on z just because of its different constant DE state parameter w(z). No 
constraint holds on H(z), any more; thus, we can take the same H value in all models at z = 0. 

Notice that the w(z) dependence so defined is quite different from w(z), the DE state equation 
in A. This outlines a significant experimental risk, as observers fitting data against constant w 
models, would find w(z), instead of the actual DE state equation w(z). 

Here we extend the test in [7] to a set of models, sampling cosmologies consistent with 
data [U [9] , and confirm that the weak requirement works up to the required k (see also [7J ) . 
We also test the requirement including hydrodynamics, confirming spectral regularities up to 
k ~ 10/iMpc -1 (see also |10|). Further hydro dynamical tests are in progress. 

A side result we shall outline concerns the rate of star formation. The size of the box 
used allows to follow such process in a rather approximate way. We however tested that our 
simulations produce an amount of star which is not unreasonable. Previous tests in boxes of 
comparable side admittedly found more substantial difficulties [11]. The point is then that, in 
spite of the proximity of the cosmologies treated, we find a significantly different star production. 
If this property is confirmed by simulations in boxes apt to follow in detail galaxy formation, it 
could envisage a new important feature to trace DE nature observationally. 

2. N— body tests for cosmologies consistent with current data 

Using available data, the WMAP team [51 [9] tried to constrain the coefficients wq and w a in 
a polynomial expression w(a) = wq + (1 — a) w a for the DE state parameter. However, the 



likelihood ellipses (Figure [T]) have significantly shifted from WMAP5 to WMAP7. The setting 




Figure 1. Likelihood ellipses on the w Q -w a plane from WMAP5 and WMAP7 data releases. 
Black triangles indicate the A models considered. State equations beyond the w Q = —w a line 
should be modified at high-z. The dotted lines cross on the WMAP7 best fit. 
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Figure 2. Each box refers here to a single A model in Figure [T] (black triangle) and related 
W(z) models (color polygones). Notice that, in most cases, distances between colored polygones 
are smaller than their distance from the black triangle. This outlines the possibility of a serious 
bias in data analysis, if data are tested by assuming constant-u> cosmologies. 



of the six A models considered is shown in the Figure. The other parameters are consistent with 
both WMAP5 and WMAP7 data: Q m = 0.274, h = 0.7, a 8 = 0.81 and n s = 0.96. 

We test the spectra of A models vs. the corresponding W(z). All simulations start at z = 24, 
from realizations fixed by an identical random seed. We use the pkdgrav code [6], modified 
to deal with any variable w(a), as in [7J. In a box of side Lb ox = 256/i~ 1 Mpc we set = 256 3 
particles and use a gravitational softening e = 25/i _1 kpc. We have 4 auxiliary models W(z) for 



z = 0, 0.5, 1, 2 (see Fig. [2]), for each A model, run just down to the redshift where they are 
tested. For all models, mass functions agree with Sheth & Tormen [12] predictions. 

The first A cosmology considered, consistent with WMAP5, is outside the 2— a curve for 
WMAP7. The significant shift of the ellipses is not due to the fresh CMB inputs. Rather, 
omitting the distance prior, as well as the whole SDSS data [HJE], surely had an impact on it. 

Figure [2] shows that the distance between W models, on the w a = line by definition, is mostly 
smaller than their distance from A. Distances depend on w Q and are smaller for the central w Q 
values. There is therefore a possible observational problem, as outlined in the Introduction. 

Model spectra are worked out by FFT-ing the matter density field, computed on a regular 
grid Nq x Nq x N G (with N G = 2048) via a Cloud in Cell algorithm. 




Figure 3. Spectral discrepancies at the redshift indicated by line colors. Each box refers to an A 
model with the same ordering as in Fig.(2j Top discrepancies at z = 0; only there they approach 
2 % for some model with w a <C — 1 and for k ~ 4-5 h Mpc -1 (here however gas dynamics should 
not be ignored). Within ~ 2-3/iMpc~ 1 , discrepancies keep < 1% . 



Figure[3]then shows spectral discrepancies. Within k = 3/iMpc -1 the top discrepancy occurs 
for the two models with lowest w Q (< —1), at z = 0. Other models exhibit a nicer behavior, 
allowing to guess that the greatest discrepancies, between A and W(z = 0) evolutionary histories, 
occur for models with w Q <C — 1 . As expected, discrepancies decrease at higher z, where non- 
linearities had no time to affect low fc's. At z = 0.5 they are in the permil range, for all models. 



3. Hydrodynamical tests 

We then compare A and W(z) models, also including baryon physics. Hydro simulations are 
much more time demanding and here we present results for one model: w Q = —0.8, w a = —0.754 ; 
top right in Figures [2] and [3j The assigned Q m value comes from Qf, = 0.046, il c = 0.228. 

The program used is gasoline, a multi-stepping, parallel Tree-smoothed particle hydro 
(SPH) n-body code [13]. It includes radiative and compton cooling. The star formation, based 
on [2], allows gas particles above constant density and temperature thresholds in convergent 




Figure 4. L.h.s.: Ratio between stellar and total mass in haloes. Empty and solid squares show 
the results for the 256 box; crosses for the 64 box. Results are compared with data in [18]; the 
40% scatter yields ~ ±2<r. R.h.s.: Star formation evolution in ACDM and dDE in the 256 box; 
the dashed (black) line are data fits [19], assuming a Salpeter IMF; curves are normalized to the 
z = value. The lower panel shows the ratio between star productions in dDE and ACDM . 



flows to spawn star particles at a rate proportional to the local dynamical time [15]. The program 
includes SN feedback [15J and a UV background, following [16] . We applied gasoline the same 
changes made in pkdgrav to handle dDE. Further details can be found in [7] 117]. 

All simulations use (2x)256 3 particles and most are run in a box of 256/i~ 1 Mpc. CDM (gas) 
particles have mass m c h/M Q = 6.33 x 10 10 (m b h/M Q = 1.28 x 10 10 ). The force resolution e 
(softening) is 1/40 of the intra-particle separation. For 256/i~ 1 Mpc this yields e ~ 25 /i _1 kpc 
(wavenumber k = 2ir/e ~ 150/iMpc -1 ). All simulations start at z = 24. The A model is 
selected so that W(z = 0) is ACDM. This ACDM was also run in a box of 64 /i _1 Mpc, still with 
(2x)256 3 particles. All simulation parameters scale accordingly. 

Figure H] (l.h.s.) shows the ratios M star l^haio vs. Mhaio in both 64 and 256/i _1 Mpc boxes, 
compared with predictions in [18]. The Figure shows that star formation generally follows the 
halo occupation distribution trend. Slightly too many stars form in haloes in the top Mh a i bin, 
due to overcooling. This might be fixed by including AGN effects. 

The r.h.s. of Fig. [5] shows star formation histories in 256ft, _1 Mpc boxes. Data [19] are also 
shown, with ±20 % error bars, approximating 2-cr's. Our outputs are shifted upwards by a factor 
1.5, so to account for an overall deficiency in star formation. Furthermore, observations and 
simulations peak for z ~ 2.3 and ~ 1.5 : stars form later in simulations. These discrepancies 
arise because stars do not form in halos with < 200 gas particles [20J, corresponding to a mass 
~ 1.3 x 10 12 h~ 1 M@ in our simulations. In spite of that, at our resolution level star formation 
seems described well enough to allow spectral comparison on the scales relevant to us. 

The bottom frame shows the ratio between star formation in A and ACDM. The decrease by 
~ 20% arises from a change of w(z), yelding a different evolutionary history, within an identical 
conformal age. Future work will seek the cause of this unexpectedly large shift. 

Let us now come to spectral comparisons between A and W(z) models at z = and 0.5, 1, 2 . 
They are shown in Figures [5j For the model considered, regularities persist when gas dynamics 
is important, although residual discrepancies (mostly of the order of a few permils) are greater 
than in n-body simulations. The largest discrepancy is found in the gas spectra, even if the gas 
spectrum distortion at z = 0, attaining ~ 5 %, apparently does not imply a significant distortion 
in the global spectrum. 

Altogether we conclude that even the inclusion of hydrodynamics keeps the discrepancies 
between auxiliary and true model in the permil range, for the model considered. 
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Figure 5. Comparison between A and W(z) spectra, for n-body (DMO) and hydro (DMG) 
simulations (various components shown). Top panels are for z = (l.h.s.) and z = 0.5 (r.h.s.); 
bottom panels for z = 1 and z = 2. 
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